Charge asymmetries in 77 



£ + £ + v's (£ — fj,, e) with polarized photons in the 
Standard Model. 
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Abstract 

It is shown that in reaction 77 — > £ + £~ + u's at \/s > 200 GeV with polarized photons, large and well observable differences 
arise in the distribution of positive and negative charged leptons (£ = fi^ , e^), (charge asymmetry). The modification due to 
the contribution of the cascade processes with intermediate r-lepton in 77 — > W + u's reaction is taken into account. This 
charge asymmetry is potentially sensitive to effects of physics beyond the standard model at the anticipated luminosity of the 
Photon Collider mode of the future international linear collider. 



I. INTRODUCTION 



The Photon Collider (PC) option of the planned Inter- 
national Linear Collider (ILC) (see e.g. [1], [2]) will offer 
a specific window for the study of new effects in both 
Standard Model (SM) and New Physics. In particular, 
it is expected that the charge asymmetry of leptons, pro- 
duced in the collision of neutral but highly polarized col- 
liding particles 77 — > £ + £~ + neutrals (where I = fx, e), 
can be a good tool for the discovery of New Physics ef- 
fects. With this aim the study of such asymmetry in SM 
is a necessary step for both better knowledge of SM and 
understanding of background for New Physics effects. 

In this paper we study the SM process, in which neu- 
trals are neutrinos and the main (but not single) mech- 
anism for charged lepton production is given by 77 — > 
W + W~ process with subsequent lepton decay of W. The 
latter process, (7(77 — > WW)Br(W — > (iv) = 8.8 pb, will 
ensure very high event rate at the anticipated integrated 
luminosity of ILC (100 fb^ 1 ), about 10 6 events per year. 
The charge asymmetry here appears due to transforma- 
tion of initial photon hclicity into distribution of final 
leptons via P-violating but CP-preserving leptonic decay 
of W. In the following we consider the particular case 
I = fi for definitcness. The considered effects are identi- 
cal for electrons and muons. So that, absolutely the same 
asymmetry will be observed in e + e~ , e + fj,~ , fi + eT dis- 
tributions. All these contributions should be added for 
a complete analysis. This will enhance the value of the 
cross section for 77 — > + v's from 1.2 to 4.8 pb. 

In the main body of the paper we consider the col- 
lision of a photon with helicity Ai moving in the pos- 
itive direction of the z axis with a photon of helicity 
A2 moving in the opposite direction. This initial state 
is denoted as 7a x 7a 2 with Ai = ± (left or right cir- 
cular polarization). For example, the initial state with 
Ai = +1, A2 = —1 is written as 7+7-. With this choice 
of the positive direction of the z axis we define the lon- 
gitudinal momentum p\\ = p z and the transverse mo- 



mentum pj_ = yp 2 +Py- For definiteness, we present 

most of the results for monochromatic photon beams at 
Vs 77 = 500 GeV (Ej = 250 GeV). The above defini- 
tions will be slightly modified when discussing the effects 
due to the non-monochromaticity of photon beams in the 
future Photon Collider. We start our numerical calcula- 
tions with the CompHEP package [3] and then switch to 
the CalcHEP package [4] which allows one to take into 
account the circular polarization of the initial photons 
and choose different random seed numbers for the Monte 
Carlo (MC) generator which is necessary for an estimate 
of the statistical inaccuracy of future experiments. 

The observable final state with W + fi or two muons 
with missing transverse momentum carried away by neu- 
trinos can appear either via processes 



77 -> Wfiv (77 -> v^v^) , 

or via cascade processes like: 

77 — > W + U T T~ 



= W + V T I1 V^Vt 



(1) 



(2) 



77 — > T V T V T T T 

i_ I 

11 V^Vr [l+V^Vr 



(J, ll + V^V T V^U T V T V T , (3) 



in which six or eight particles are present in the final 
state. 

To reduce CPU time in the Monte Carlo event gen- 
eration, 10 6 events for each channel, we obtain the es- 
sential part of the results for the 77 — » M /± /i =F + v's 
process, not taking into account issues related to the 
reconstruction of the W . The analysis of this process 
allows us to extract the main features of the effect of in- 
terest, i.e. the difference in the distributions of [i + and 
/i~ at fixed photon helicities (global charge asymmetry). 
We show that the additional diagrams that contribute to 
77 — > £ + £~ + neutrals give negligible contribution to the 
cross section and charge asymmetry. 
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We start with the description in Section II of the gen- 
eral features of the effect, neglecting cascade processes. 
We classify the diagrams contributing to 77 — > W \xv 
and 77 — > fjr pTv^v^ according to the different topolo- 
gies and give an approximate analytical estimate of their 
relative impact on the cross section (Sec. II A). These 
estimates allow us to present qualitative explanation of 
the appearance of the charge asymmetry (Sec. II B). 

We then introduce suitable variables for the description 
of the global asymmetry, i.e. the difference in distribu- 
tions of n + and fi~~ in the processes 77 — * W ± p ± + v's 
or 77 — > n + + v's (Sec. II D). In this very section 
we describe the cuts applied to the observed particles. 
In Sec. HE we discuss a computational method used to 
estimate a lower bound on the statistical uncertainty of 
future experiments. This estimate is obtained directly 
by the repeated Monte Carlo simulations with an antici- 
pated number of events. 

Section III is devoted to detailed description of the 
global charge asymmetry of leptons in the process 77 — > 
Win; in monochromatic 77 collisions. 

The accurate calculation of cascade processes with six 
or more particles in the final state is a computationally 
challenging task with available software. Since we use 
CompHEP/CalcHEP packages which don't fix the he- 
licity of final states while the discussed effects strongly 
depend on the helicity, the direct use of existing software 
for tau decay simulation like TAUOLA [5] is not possible 
here. In Sec. IV we construct reasonable approximations 
in the description of cascade processes (3). The detailed 
analysis of the modification of momentum distributions 
allows to find that the inaccuracy implemented by the 
mentioned approximation in the final result is within the 
estimated statistical uncertainty of future experiments. 

In Sec. V we discuss the total observable asymmetries. 

High energy photons will be produced at the Photon 
Collider through Compton back-scattering of laser pho- 
tons from high energy electron or (and) positron beams: 
the photons will not be monochromatic but will demon- 
strate an energy and polarization distribution. The high 
energy part of this spectrum will mainly include photons 
with definite helicity A.; close to ±1 [2]. We analyse the 
influence of initial photon non-monochromaticity on re- 
sults in Sec. VI. 

The correlative asymmetry in fj, + and fj," momenta in 
each event of 77 — > \i + [i~ + v' s is expected to provide 
more information in the search for effects of physics be- 
yond the SM. We discuss it in Sec. VII. 

We conclude and summarize the obtained results in 
Sec. VIII. 

Preliminary (and incomplete) parts of this work were 
reported earlier [6]. 



II. GENERAL FEATURES 

The SM cross section of 77 — > W + W~ at center of 
mass energy greater than 200 GeV remains almost con- 




FIG. 1: Classes of tree level Feynman diagrams contributing 
to 77 — > W + fi~i>, (a-c). The grey blob in (a) represents 
diagrams with W exchange with trilinear yWW coupling and 
the diagram with quartic jjWW coupling. 



stant at the asymptotic value a ~ 8ira 2 /M^ ~ 80 pb 
and practically independent on photon polarization [8], 
see the formulas in Subsection II B. At y/s > 200 GeV 
this cross section is more than ten times larger than 
the cross section of W production in e + e~ mode. It 
will ensure very high event rate at the anticipated lu- 
minosity. The distributions of W + and W~ bosons in 
77 — > W + W~ process are identical (charge symmetrical 
distribution), their polarizations are determined by the 
polarization of initial photons. The distribution of muons 
in subsequent decay of polarized is asymmetrical due 
to P non-conservation with CP conservation in the SM. 



A. Diagrams 

In this section we classify all tree level diagrams de- 
scribing the process 77 —> W^p^v and 77 — > /i + pTvv 
in classes according to their topology (a similar classifi- 
cation was given also in Ref. [7]). For each topology we 
give an analytical estimate of its asymptotic contribution 
to the total cross section at s ^> M$y, identifying in each 
group the 2 — ► 2 dominant subprocess and assuming for 
the SM gauge couplings g 2 ~ g' 2 ~ e 2 = Ana. The nu- 
merical Monte Carlo results, supporting these estimates, 
are presented in the next Sections. 

The processes 77 — * W^p^v are described by seven di- 
agrams, which we divide in three classes, shown in Fig. 1: 

(a) Three double-resonant diagrams (DRD) of 
Fig. 1(a) describe WW pair production with 
subsequent decay. Their contribution to the total 
cross section is ~ a 77 ^wwBr(W — > pv) ~ 
{a 2 /M^)Br(W^pu). 

(b) Two single-resonant diagrams (SRDW) of Fig. 1(b) 
with W exchange in ^-channel contribute to the 
total cross section a s ~ aa^^ww ~ ( a /^w)- 
The relation between this contribution and DRD 
contribution is o~ s /ad ~ a/Br(W — > \xv). 

(c) Two single resonant diagrams (SRD/x) with 
lepton exchange in i-channel (gauge boson 
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FIG. 2: Additional tree level Feynman diagrams contributing 
to 77 — > [i + [i~vv: the grey blob in (d) represents diagrams 
with fifj, fusion to Z and diagrams with Z radiated by an 
external [i line. 



decay and permutations of external fermion lines (with 
the same estimates as above) and two additional types of 
diagrams shown in Fig. 2. 

(d) Six diagrams with radiation of Z boson in the 
process 77 — > fi + fi~, Fig. 1(d). The asymptotic 
contribution is az ~ cta 11 ^ ll ^aBr(Z — > vv) ~ 
(a 3 /s)Br(Z -> vv). 

(e) Two multi-peripheral non-resonant diagrams 
Fig. 1(e) with a n ~ a 4 /M w . 



bremsstrahlung) , Fig. 1(c). The contribution to 



the total cross section is a 



'77— i/MfJ, 



(ays), 



therefore a stl ja^ ~ [a/Br(W -> fiv)](M^/s) 

The process 77 — > ji/ivv is described by the diagrams 
Fig. 1 with the addition of lines describing the W ^> fiv 



J 



B. Qualitative picture 

The above analysis shows that the bulk of the cross 
section is given by the diagrams containing the process 
77 — > W + W~ with subsequent decay of W bosons to 
leptons Fig. 1(a), DRD diagrams. 



Denoting by p± W's transverse momentum, the 77 — * W + W differential cross section can be written as 

8na 2 



da = dcr n P + X 1 X 2 dT a , 

da n P = a w 
dr a = aw 



aw = 
(16 



4M^ 



M w ' 
■3x 2 )M w 



dt 



dpi 



(3 + 8x)x 



s 

3x 2 1 



[32(pl + M w y 



xM w 



32(p 2 



M w ) 



+ 



(3 + 8x)x 



2{p\+M w y 32(pi + M 



6AM W 
3x 2 



dt, 
dt, 



with total cross section (v = yl % ) 
a = awv ■ < 1 



1 (1 ) log — 

16 16 V 2 ; 16i; 5 1 



xAiA 2 
16 



-19 



(8 - 5x) , 1 
log — 



(4) 



(5) 



From the analysis of these equations we can see that 
at yfs > 200 GeV the 77 -> W+W~ differential cross 
section practically does not depend on photon polar- 
izations, as for standard QED process, total cross sec- 
tion is practically energy independent, at the plateau 
value aw — 8ira 2 /M w . Moreover, the W's are pro- 
duced mainly in the forward and backward directions, 
with average transverse momentum ~ Mw (distribution 
<xl/(pl+M w ) 2 ). 

As shown in Ref. [9], in this process we have an ap- 
proximate helicity conservation. For p± — 0, the he- 
licity of W ± moving in the positive direction of z axis 
is Xwi = Ai, irrespective of the charge of the W; the 
same holds for the W moving in the opposite direction: 
Xw 2 — M- These identities do not hold for p± ^ and 
become less and less accurate with increasing values of 
p± . Since in our process the cross section is concentrated 
at small values of p±, we have an approximate helicity 
conservation: Xwt ~ Ai, and Xw 2 ~ A2, both for W + 



and W~. 

Now we can qualitatively understand the origin of 
charge asymmetries. Let z'-axis be directed along W 
three-momentum and e « Mw/2 and p z i be the energy 
and the longitudinal momentum of \i in the W rest frame. 
It is easy to calculate that the distribution of muons from 
the decay of W with charge e = ±1 and helicity A = ±1 
in its rest frame is cx (e—eXp z i) 2 (the transverse momenta 
of muons are distributed roughly isotropically relative to 
W momentum within the interval p± < mw/2). In other 
words, the distribution of muons from decay has a 
peak along W momentum if the e • Xw = — 1 and oppo- 
site to W momentum if e • Xw — +1- These distributions 
are boosted to the distributions in the 77 collision frame. 
For example, for a collision of photons in a 7_7- initial 
state, the fj," are distributed around the upper value of 
their longitudinal momentum (in forward and backward 
direction), while the // + are concentrated near the zero 
value of their longitudinal momentum. At the same time, 
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this boost makes the distribution in p± wider in the first 
case and narrower in the second case. 



C. Cuts 

It is natural to expect that the relative size of New 
Physics effects will be enhanced with the growth of trans- 
verse momenta of observed particles. This is the main 
reason why we study the dependence of observed effects 
on the cut in p± . Namely, we impose cuts on the trans- 
verse momenta of observed charged particles, pj_ and on 
the scattering angle 



P± > Pin > 
6 <9<ir-6 , 9 Q =pl/2E. 



(6) 



The cut (6) is applied to each observed particle and to the 
total transverse momentum for the sum of momenta of all 
observed particles; the cut for escape angle is applied to 
all the observed particles. We consider the dependence of 
all studied quantities on the pj_ up to p°. = 140 GeV. 
This cut also mimics limitations from the detector in the 
future experiment. 

These simultaneous cuts allow to eliminate many back- 
grounds (since charged particle (s) with missing trans- 
verse momentum greater than p± should have the es- 
cape angle greater than 1p\J\[s > a ). In particular, 
all pure QED and QCD processes are eliminated by these 
cuts, since they cannot provide large missing transverse 
momentum. 

If it is not otherwise specified, in the following we set 
Pj_ = 10 GcV, #o = 20 mrad, and monochromatic pho- 
ton beams at = 500 GeV {E 1 = 250 GeV). 

The effect of cuts with non-monochromatic photons is 
studied in sect. VI. 



D. Variables for description of global asymmetry 

The global asymmetry variables are described by the 
difference in distributions of ^+ and in the processes 
77 — > W + fi~ + f's and 77 — > W~ n + + v's (or in 77 — > 
H + fi~ + i/'s). For definiteness, we calculate all quantities 
only for the case when negatively charged particle (W~ or 

) is in the forward hemisphere (p\\ > 0). In the study 
of the dependence on p\^ , we will label all the quantities 
by the argument (Pj_„)- 

A suitable measure for the longitudinal (At) and trans- 
verse (Ay) charge asymmetries are the relative differ- 
ences of corresponding momenta distributions for nega- 
tive and positive muons: 



J pn da — J ptda 



Ip\\ da + Jp+da 



A T 



j p ± da — J p~\_da 
J p^da + J p\da 



(7) 



It is useful to define also mean values of longitudinal pj 
and transverse p"J_ momenta of or ^+ 



I ptda 



E. 



jmax 



J da' 



JTrr 



Jp^da 



E, 



■ymax 



J da- 



(These definitions are written in the form which is useful 
for non-monochromatic case as well.) 
Due to CP symmetry of the SM, 

d°—(Pn+,Pn-) =da++(p ll -,p t , + ) 
da+-{p^+,p^-)^da^ + (p^-,p^). (9) 

(Here subscripts + and — at cross section label initial 
photon helicities). 1 . In particular, the distributions of 
/i~ and [i + in the forward hemisphere for 7+7- collision 
reproduce the distributions of /i + and /i~ in the backward 
hemisphere. Therefore, in all the cases the asymmetries 
Al (determined by Eq. (7)) change signs in each hemi- 
sphere when the helicity changes to opposite. For the 
7-7+ collisions these asymmetries in forward and back- 
ward hemispheres have opposite signs. These symmetries 
break if any CP- violating interaction is present. 

Total cross sections of the processes 77 — > VF + /i~ + v's 
and 77 — ► W~ fi + + v } s coincide at each initial photon po- 
larization. However, in accordance with above discussed 
qualitative picture, applied cuts reduce these cross sec- 
tion in different way. So that, it is useful to define relative 
value of this difference in dependence on cut variable, 



Aa( P i M ) 



(fda(W-»+)-fda(W+»-)) {p±>p%) 
(jda(W-^)+jda(W+^)) {p±>p ^ 



) 

(11) 

and the fraction of the total cross section left by the cut 
in pl^. 



a ± (p^) = j da(W T ^)\ 



(p±>p1J- 



(12) 



E. Estimate of statistical uncertainties 

MC calculations simulate an experiment and have 
some statistical uncertainty Smc- This uncertainty value 



1 One can write the differential distribution in the reaction 77 
+ v's as 



da 



A + B\\ +CX 2 +-DA1A2 



(10) 



That is another form of eq. (9) with J BcPp^+cPp^- = 0, 
J C 'd 3 'p M + d 3 'p fl - = 0. The weak dependence of the cross sec- 
tion on the photon polarization means that / Dd 3 p fi+ d 3 p^- < 
f Ad 3 p^+d 3 p^- . Our subsequent analysis based on momentum 
distributions shows that on average |U|~|S|~|S|~|A|. In 
the following paragraphs and sections we will not make use of 
the form in Eq. (10). 
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for the integral characteristics like (7) cannot be pre- 
dicted simply from general reasons. To find this uncer- 
tainty we repeated our MC calculation with anticipated 
10 6 number of events five times for different random num- 
ber inputs for MC generator. Additionally we consider 
as an independent input the set of observations obtained 
by simultaneous change Ai, A2 — ► — Ai, — A2, p~ <-» n + 
(this change should not change distributions due to CP 
conservation in SM), in whole it corresponds ten repeti- 
tions of "MC experiment" in a sum. These sets of data 
were an input for standard generation of Monte Carlo 
inaccuracies Smc- 

Since the adaptive MC is used with CalcHEP, it is nat- 
ural to expect that the statistical uncertainty of future 
real experiment 5 S J^ > Smc- Therefore, below we omit 
subscript MC, having in mind that our numbers give es- 
timate for statistical uncertainty from below. 



III. GLOBAL ASYMMETRIES IN THE MAIN 
PROCESS 77 -f W ± ^v. MONOCHROMATIC 
CASE 



We start with the study of asymmetry neglecting the 
cascade channel and supposing photon beams monochro- 
matic and completely polarized. First, we present the 
distributions d 2 a/ (dp\\dp±) of muons in the (pu, p±_) 
plane, at different photon polarizations, in Figs. 3. These 
figures show explicitly strong difference in the distribu- 
tions of negative and positive muons as well as strong 
dependence of distributions on photon polarizations. 
Therefore, the charge asymmetry in the process is a 
strong effect. 

Table I presents obtained average momenta for the neg- 
ative and positive muons and corresponding asymmetry 
quantities (7) together with their statistical uncertainties 




FIG. 3: Muon distribution in 7-7- — > Wfiv (upper plots) 
and in 7+7- — > Wfu> (lower plots), left column fj,~ , right 
column n + . 
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FIG. 4: Dependence of cross sections on cut p c ± ^ for main 
process. 

(in percents) for different cuts Pj_ . One can see that the 
values of asymmetry are typically 20-50%. 

We have checked that with the change of sign of both 
photon helicities mean muon momenta for negative and 
positive muons change their places (within statistical ac- 
curacy) so that the quantities A^t change their signs 
with this change of polarization. 

The longitudinal scale of distributions in momenta is 
determined by initial photon energy while the transverse 
scale is determined by the W mass. Hereupon mean 
transverse momenta are usually smaller than longitudi- 
nal. 

Besides, for the collision of photons with identical he- 
licity at the growth of cut p±^ the cross sections for pro- 
duction of positive and negative muons become different, 
due to discussed charge asymmetry (remind, we discuss 
only events with negative muons or Ws flying in for- 
ward hemisphere). For the collision of photons with op- 
posite helicity these cross sections should coincide since, 
for example, for 7-7+ collision forward hemisphere for 
/i~ production realize absolutely the same distribution 
as backward hemisphere for /z + production. This effect 
is clearly seen from Figure 4, where we present the p c ±li 




20 40 60 80 100 120 140 100180 200 220 240 




P^' GcV 



FIG. 5: Relative cross section differences in dependence on 
cut p^ M value. 
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10 


7-7- 


0.606 


0.29% 


0.201 


0.55 % 


+0.501 


0.57% 


0.333 


0.61% 


0.159 


0.28% 


+0.355 


0.44% 


7+7- 


0.223 


0.74% 


0.609 


0.19% 


-0.463 


0.47% 


0.164 


0.08% 


0.262, 


0.31% 


-0.231 


0.76% 


40 


7-7- 


0.593 


0.39% 


0.273 


0.20 % 


+0.370 


0.47% 


0.378 


0.64% 


0.241 


0.62% 


+0.222 


1.07% 


7+7- 


0.296 


0.64% 


0.637 


0.25 % 


-0.366 


0.66% 


0.239 


0.28% 


0.319 


0.25% 


-0.143 


0.31% 


140 


7-7- 


0.402 


0.68% 


0.242 


0.14% 


+0.249 


1.45% 


0.697 


0.11% 


0.621 


0.04% 


+0.057 


0.95% 


7+7- 


0.253 


0.81% 


0.489 


0.27 % 


-0.318 


1.33% 


0.672 


0.12% 


0.660 


0.05% 


-0.009 


5.75% 



TABLE I: Charge asymmetry quantities and their statistical uncertainties for the process 7a x 7a 2 — * Wfw. 



dependence for cross sections (12) at different initial pho- 
ton polarizations. 

The corresponding cross section differences, Eq. (11), 
are shown in Figure 5 without averaging over realizations. 
The difference in cross sections for 7-7- and 7+7+ colli- 
sions is due to our choice of events with negative particles 
in forward hemisphere. In this case the cross sections 
coincide at small Pj_ , while at high p c , one of them 
becomes larger and larger in comparison with the other, 
thus Ac goes to one (one more demonstration of trans- 
verse asymmetry). For the case with opposite photon 
helicities the deviation of Act from zero shows high sta- 
tistical uncertainty, due to the very low value of the cross 
sections (low counting rates) at p±^ > 120 GeV. Similar 
dependencies with roughly the same characteristic values 
of p\ remain valid even at higher collision energies (for 
example, at ^fs = 2 TeV) since the scale of this depen- 
dence is determined by the W boson mass and not by 
the total energy. 

IV. GLOBAL ASYMMETRIES IN CASCADE 
PROCESS WITH INTERMEDIATE r 

The observable final state with two muons or W + \i 
with missing transverse momentum carried away by neu- 
trinos can appear either via processes 77 — > fj, + fj,~ i/^D^ 
(77 — > Wfiu) or via cascade processes with r produc- 
tion and subsequent r decay (r — ► [iv ll v T ), see Eq.(3). 
The latter process enhances the total event rate (with- 
out cuts) by a value given by factor B = Br(r — + \xvv) — 
17 % for the 77 — ► Wfi+v's. Similar event rate enhance- 
ment in the process 77 — ► + v's is 2B + B 2 w 37 %. 

The accurate calculation of processes with 
six or more particles in the final state like 
77 — ► /i + pT ~ v^v^v-tV-t is a computationally challenging 
task with available software. Since r is very narrow 
particle, the diagrams without r-pole in s-channel can be 
neglected with very high precision, (<~ T T /m T ). There- 
fore, one can in principle use the results for 77 — > Wtv 
(77 — > r+ 'li'Vypp, etc.) and convolute them with the 
distribution of fi from r decay. However, the latter 
distribution depends on r polarization which cannot 
be determined definitely with CompHEP/CalcHEP in 
general case. (Generally, all helicity amplitudes for r 
production are nonzero, and for convolution one must 



consider not only diagonal helicity states but also their 
interference). 

Fortunately, the cascade process provides only a small 
fraction of the total cross section and the main contri- 
bution to the total cross section is given by the double 
resonant (DRD) diagrams of Fig. 1(a). That is the rea- 
son why in the description of the cascade contribution 
only these diagrams can be taken into account — DRD 
approximation. In this approximation the r helicity is 
precisely determined in each MC event. In sect. VA we 
will show that the inaccuracy introduced in the total re- 
sult using the DRD approximation for the cascade con- 
tribution is within the estimated statistical uncertainty, 
found for the main process. 

Note that the distributions obtained in sect. Ill de- 
scribe with high accuracy also r-distributions in 77 — ► 
Wtv processes etc. 

In the DRD approximation each r is produced only via 
W-decay, and its polarization in the rest frame of W is 
given by the SM vertex, t + W~7^(1 — 7 5 )i+ + h.c Due 
to 7^(1 — 7 5 ) factor, r helicity is opposite to that of v T , 
it is positive for r + and negative for t~ (with accuracy 
to m T /Mw) and independent on W polarization. 

For each generated event momenta of all particles are 
known. The spin vector of r is expressed easily via mo- 
menta of t and v r , p T and p v respectively, as 

±s/2, where s = ( * - f^) {+ 

\m T (p T p v )) [ - for t . 

(13) 

Denoting the momentum of // by k, the distribution of 
muons in T-decay with momentum p T and spin ±s can 
be written, neglecting muon mass, as 

/ = — g [(3m^ — Ap T k)p T k + ks ■ m T (Ap T k — m 2 T )\ aT, 

(14) 

where dT is a phase space element boosted to the lab 
frame. In the r rest frame dT = 9 (m T /2 — k) d 3 k/E^. 
Note that the sign of helicity before s in Eq. (13) disap- 
pears in the result. 

Let us discuss now the qualitative features of the muon 
spectrum given by the convolution of the r spectrum with 
the distribution in Eq. (14). One can consider r — > [ivv 
decay as a two body decay: a massless muon and the 
di- neutrino with invariant mass m vv . At given m vv , the 
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energy and 3-momentum of the muon in the r rest frame 
are £ M = p° = Am T /2 with A = 1 — (m^ u /m T ) 2 . In 
the laboratory frame where the 3-momentum of the r, 
p T is under some angle 9 relative to muon momentum 
in t rest frame, the muon momentum is evidently p^ = 
Ap T (l + cos0)/2 plus small corrections negligible in our 
discussion. Therefore, the muon distribution repeats in 
some sense that of the r but with a factor A(l + cos#)/2, 
which is usually much lower than one. In other words, 
the distribution of muons in the cascade process is similar 
in the main features to that of the r but it is strongly 
contracted to the origin of the coordinates. 

• It is useful to describe the inaccuracy of the DRD 
approximation in the description of the 77 — > Wtv cross 
section itself, 8J^ D , in dependence on the cut p\ T . The 
estimates in Sec. II A show that at the considered en- 
ergies all contributions to the cross section are small in 
comparison with that of the DRD except for the SRDW 
contribution. The interference term Re(A* DRD AsRD\v) 
is roughly of the same order of magnitude as |As_r,cw| 2 
since the DRD is large only in regions of the final phase 
space corresponding to the W resonances, while the other 
contributions do not have these peaks. The numerical 
value of this inaccuracy is obtained by direct comparison 
of this cross section, calculated with all diagrams, and 
those for DRD diagrams with MC simulation. 

At p\ T = 10 GeV wc find that for the 77 — > Wtv 
process the SRDW contribution itself is about 5 % of 
the DRD one, and the interference of this contribution 
with DRD is destructive so that the DRD contribution 
differs from total cross section only by about 1%. This 
difference naturally grows with increasing values of p\ T . 

More important for us is the inaccuracy of DRD ap- 
proximation in the description of asymmetry quantities 
(7). Table II presents value of inaccuracy Sjj RD (p° Lt ) 
introduced by DRD approximation in the description of 
77 — ► Wtv process at different cut values of r transverse 
momenta p\ T . 




FIG. 6: Muon distribution in cascade processes 7-7- — > 
Wfivvv (upper plot) and 7+7- — > W jivvv (lower plot), left 
- £t~, right - M + 



pi T> (GeV) 
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10 


0.9 


2.3 


0.7 


3.45 


40 


1.5 


3.6 


2.1 


4.1 


80 


1.9 


5.6 


4.2 


7.7 


120 


5.7 


5.3 


4.4 


31 



TABLE II: Inaccuracy of DRD approximation 5dhd(p1t) f° r 
Al,t at different p c ±T for r production. 




FIG. 7: Total muon distribution in 7-7- — > Wfx + v's (upper 
plot) and 7+7- — > Wfi + v's (lower plot); left - u~ , right - /i + 



Large values of the relative quantity 6 T at large pj_ T 

appear in the case when the absolute value of Aj!~ ^ is 
negligibly small. One can see that this inaccuracy grows 
with increasing the value of the cut, see Fig. 5. However 
with this increase also the fraction of cascade process 
within the total process becomes smaller and smaller (see 
discussion at the end of next section) . 

For the processes 77 — ► fi^T^vv and 77 — > 
t + t~vv the inaccuracies of DRD approximation 

are sT drd(p1t) = s rmD(P±r) and 5 drd(p1t) = 
2^(P1 T )- 

V. TOTAL ASYMMETRIES 

The resulting distributions include the complete tree- 
level results of 77 — ► W [iv and DRD approximation for 
cascade contribution. 

Figures 7 show the total observable distributions of 
muons, i.e. the sum of distributions of muons in 77 — > 
W [iv and 77 — > Wtv — > Wfivvv, and Table III presents 
the corresponding total asymmetry quantities for p\^ = 
10 GeV. Comparison with Fig. 3 and Table I shows that 
the cascade process introduces a change in the shape of 
muons distribution only at small momenta and its con- 
tribution reduces the asymmetry parameters A-l.t in av- 
erage by about 3 % only. 

The plots in Fig. 8 show the dependence of asymme- 
tries Al and At on p\^- The longitudinal charge asym- 
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-- polarization 
-+ polarization 




150 200 
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150 200 
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FIG. 
right 



8: The p c ±ll dependence of asymmetry. Left plot - Al, 
plot - At, black lines - for 7-7-, gray lines - for 7-7+ 



metry remains large even with large cuts, while the trans- 
verse charge asymmetry diminishes with p c j_ growth. In 
particular, for 7+7- collision at pj_ > 120 GeV the 

quantities P^ and P^ practically coincide, giving neg- 
ligible At, right plot, with naturally high statistical un- 
certainty in this small quantity. 

The Pj_ dependence of resulting cross sections for dif- 
ferent photon polarizations is shown in Fig. 9. The differ- 
ence in curves for (++) and ( ) initial states arise be- 
cause of our charge asymmetric selection of events, with 
negative particles flying in the forward hemisphere (see 
Sec. II D). 



A. Inaccuracy of DRD approximation for result- 
ing asymmetries 

Let us denote by &drd{p Lu) the inaccuracy of the 
DRD approximation for the resulting asymmetries, by 
£ cosc (pj_ ) — the inaccuracy of DRD approximation for 
the description of the cascade process itself, like that 
given in the Table II, and by 



R(p[ 



"{Ply) 



(15) 



the relative contribution of cascade /x in the total cross 
section, all - in dependence on cut for muons p^. Nat- 
urally, 



Stot 



(p%) = R(pl„)6 casc (p%). 



(16) 



Due to the contraction of the distribution of muons 
produced in the cascade process in comparison with that 



7ax7a 2 


Pl 


Pl 


A L 


P~ 


P+ 

T 


A T 


7-7- 


0.548 


0.164 


+0.539 


0.311 


0.142 


+0.374 


7+7- 


0.199 


0.513 


-0.440 


0.152 


0.232 


-0.207 



of the parental r (see Sec. IV), and consequently with that 
of /i in the main process, the relative contribution of the 
cascade /i in the total cross section -R(f>j„) falls rapidly 
with the growth of the cut p5_„, as can be seen in Fig. 10. 

At p c ± = we have 5 casc (p c j_^ = (5 casc (j/j_ T ). With 
the numbers given by Table II and Fig. 10 one can see 
that the inaccuracies (16) are lower than the expected 
statistical uncertainty of future experiments (Table I, 
first two lines). With growth of pj_ the inaccuracy 
S casc (p c LtJ ), similar to that given in the Table II, in- 
creases, but the cascade contribution Fig. 10, 
decreases faster. Therefore the resulting inaccuracy in- 
troduced by the DRD approximation in Eq.(16) is well 
within the expected statistical uncertainty of future ex- 
periments, Table I, at each cut on transverse momentum 
for the 77 — > W [i + v's process (and it is within the ex- 
pected statistical uncertainty of future experiments for 
77 — ► n + n~ + v's process). 



VI. EFFECT OF 
MONOCHROMATICITY 



PHOTON 



NON- 



At the PC photons will be non-monochromatic with 
spectra peaked near the high energy limit E™ ax . More- 
over, due to the finite distance between the conversion 
point (CP) and the interaction point (IP) and also due 
to rescatterings of laser photons on electrons after the 
first collision, photon spectra are even non-factorizable. 
Fortunately in their high energy part, (E 1 > E™ ax /\/2), 
these spectra are factorizable with a high precision and 
these photons have a high degree of polarization. More- 
over, the form of the effective spectra in this region is 
described with high accuracy with the aid of only one ad- 
ditional parameter, independent from the details of the 
experimental setup, while the polarization is the same 
as for pure Compton effect [10]. The luminosity of the 
Photon Collider is normalized for this very region only. 

The low energy part of the effective photon spectrum, 
depends strongly on the details of the experimental setup 
which may change during the construction process of the 



O 'total Pb 



7-7- 

7 + 7_ and 7_7+ 
7+7+ 




TABLE III: Resulting asymmetry quantities. 



20 40 60 80 100 120 140 160 pf GeV 



FIG. 9: The p c ±ll dependence of resulting cross sections for 
different photon polarizations. 
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FIG. 10: Relative contribution of the cascade process R(p±^) 
as defined in Eq.15 vs 




scmble those presented in Fig. 3 with additional maxi- 
mum at low energies. Table IV shows the corresponding 
asymmetry quantities. These values are slightly smaller 
in comparison to the monochromatic case, but they are 
still large enough and replicate in main features the val- 
ues in Table I with approximately the same statistical 
uncertainties. 



VII. 



CORRELATIVE ASYMMETRIES IN 77 

+ v'S 



The charge asymmetry in relative distributions of pos- 
itive and negative muons in each event can be a more 
useful instrument to hunt for the New Physics (but with 
lower counting rates). A simple analogy in tcrmsOO of 
charge symmetric variables, is provided by transverse 
momentum and invariant mass distribution. The global 
asymmetry distribution corresponds to that in transverse 
momentum while the correlative asymmetry distribution 
corresponds to that in the effective p + p~ mass. The 
latter is sensitive to the existence of possible resonance 
states, which cannot be seen in global asymmetries. As 
an example we present distribution in k = p+ + p- in its 
longitudinal and transverse component, Fig. 13. 

In the case of charge symmetry this distribution would 
be centered around the point (fcii,A;_i_) = (0,0). This fig- 
ure exhibits strong effect of charge asymmetry. 

The first problem for numerical analysis here is to 
find some representative variables in 5-dimensional space 
of observables p+p~. We consider three representative 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



" natural" dimcnsionlcss variables for 77 
cess: 



[V \i vv pro- 



FIG. 11: The "realistic" photon spectra, used in our calcula- 
tions. 



ILC. Therefore, in our simulations, we used a photon 
spectrum composed of two parts as shown in Fig. 11. At 
E 1 > E™ ax /y/2 we used the approximation from Ref. [10] 
with p = 1 and x = 4.8 with polarization for ideal Comp- 
ton effect [2]. In order to imitate the low energy part 
of the spectrum (at E 7 < E™ ax /\/2) we used spectra 
from [2] for the case when the IP and CP coincide (p = 0) 
and consider these photons to be unpolarized. 

In this section we denote by 7_ an initial photon state 
obtained in collision of laser photon with helicity P c = +1 
and an electron with mean double initial helicity 2A e = 
0.85, which gives A 7 = — 1 for photons with E 1 = E™ ax . 
In this case mean polarization of photons with E~ f > 
E™ ax /y/2 is also negative but its absolute value is lower 
than one as it is described in [2]. At E 7 < E™ ax /V2 
we treat these photons as non-polarized. The state 7+ is 
defined in the same way. 

The resulting distributions of muons are presented on 
Fig. 12 for the case when incident electron energies are 
250 GeV and the laser parameter is x — 4.8. They re- 



v = 



4(p 2 



-1+ 



Pl- 



) 



4(P| + 



4 



M 2 W 

4(P||+e+ - 
w = 



M 2 W 



Ml 



(17) 




FIG. 12: The distributions of muons calculated with "realis- 
tic" spectra distribution. Upper plots - 7-7-. Lower plots - 
7+7-. Left - p~ , right - p + 
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0.10% 
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0.81% 
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0.24% 


0.338 


0.08% 


-0.321 


0.43% 


0.200 


0.09% 


0.236 


0.16% 


-0.082 


0.42% 


40 
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0.375 


0.52% 
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0.16% 
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0.65% 
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0.51% 
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0.319 
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-0.067 


0.82% 


80 


7-7- 


0.355 


0.47% 


0.208 


0.23% 


+0.263 


0.88% 


0.515 


0.08% 


0.449 


0.06% 


+0.069 


0.53% 


7+7- 


0.207 


0.38% 


0.338 


0.17% 


-0.305 


0.75% 


0.467 


0.06% 


0.483 


0.03% 


-0.017 


2.53% 



TABLE IV: Charge asymmetry quantities for "realistic" photon spectra, ^/s ee = 500 GeV. 



| A g / A K M A K | , (- ■), ^Js=500 GeV | 




-200 -100 100 200 

Klong 



FIG. 13: Distribution in fci|,fe_L, monochromatic photons, 
v/s = 500 GeV 



where e± is the energy of the /x . The asymmetry quan- 
tities are mean values of this quantities averaged over all 
events allowed by cuts. In the future study of effects of 
New Physics some other variables can be more useful. 

We present distributions in these variables in Fig. 14. 

For the 7-7- and 7+7+ collisions the distributions in 
the forward and backward hemispheres are identical. For 
these initial photon polarizations w — while the vari- 
ables u and v describe interesting asymmetries. Vice 
versa, for the 7-7+ collision the distributions in the for- 
ward and backward hemispheres can be obtained from 
each other by the exchange /i ± <-» Therefore, for 

these collisions u = v = while w describe the charge 
asymmetry, see Figure 14. 



VIII. SUMMARY AND OUTLOOK 

Let us enumerate main results obtained in this work. 

• We consider the charge asymmetry of leptons pro- 
duced together with neutrinos in the collision of polar- 
ized photons. This charge asymmetry is defined as the 
difference in the momentum distributions of the produced 
negatively and positively charged leptons, and arises be- 



cause the CP conserving weak interaction vertex makes 
the momentum distributions strongly correlated to the 
initial photon polarization. This asymmetry is observ- 
able for each fixed circular polarization of at least one 
colliding photon. 

• In particular, we present a detailed analysis of charge 
asymmetries, in the SM reactions 77 — * W £^ + v' s and 
77 — » l + Cr + v' s. with polarized photons. The method 
of observation of this effect, described in detail in the 
text, is based on a standard differential analysis of final 
state momentum distributions of the observed leptons 
with suitable applied cuts, and using well known Monte 
Carlo software for the generation of events. 

• We suggest the method for obtaining an estimate of 
the lower bound for the statistical uncertainty of future 
experiments as given by the error of the Monte Carlo sim- 
ulation at the anticipated number of events. We find that 
this uncertainty for the quantities under interest in our 



do/du I do/dv I 




FIG. 14: Distribution in u (left) and v (right) for 7-7- 
collision. Bottom: Distribution in w for 7-7+ collision. 
The 77 — > /i + fiT vv process with monochromatic photons at 
v/i = 500 GeV 
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problem is significantly larger than 1/VN (by a factor 
3 5). 

• Table I shows that the statistical uncertainty in the 
charge asymmetry is at the level of radiative corrections. 
Therefore loop corrections to the differential distributions 
and the resulting corrections to the charge asymmetries 
can be safely neglected in the analysis of this type of 
experiments to be performed at a photon collider. One 
can hope to observe the effects of radiative corrections 
only if the luminosity would be enhanced by a factor 

10 -e- 100. 

• Processes with intermediate tau lepton decays (cas- 
cade process) do also contribute to the final state with ^ ± . 
We have constructed an approximation, which describes 
cascade processes simply (based on the double resonant 
diagrams for pair production). This approximation 
describes the contribution from cascade processes to the 
observable charge asymmetries with high enough accu- 
racy, within the statistical uncertainty of future experi- 
ments. 

• Taking into account the cascade processes changes 
the charge asymmetry only weekly, the relative value of 
this contribution decreases at increasing values of the cut- 
off momentum p\ . 

• We have further shown that the non- 
monochromaticity of photons at Photon Colliders 



diminishes the considered asymmetries, but only weekly. 

• The substantial reduction of cross sections at in- 
creasing values of the cut-off momentum p c ± above 
is to be compared with the fact that, on the contrary, the 
charge asymmetries are affected only slightly by p c ± . This 
makes the charge asymmetry a very good candidate as 
optimal observable for the discovery of New Physics ef- 
fects in the processes 77 — > £ + l~ + neutrals if, as it 
is expected, the scale of the New Physics is larger than 
M w . 



Acknowledgments 

This work is supported by grants RFBR 08-02-00334-a 
and NSh-1027.2008.2. I. F. Ginzburg acknowledges sup- 
port from the Centro di Cultura Scientifica Alessandro 
Volta, Landau Network office, which allowed a visit to 
INFN Sezione di Perugia, where this work was initiated. 
This work was also partially supported, in the earlier 
stages, by the European Contract HPMF-CT-2000-0752. 
K.Kanishev is supported also by EU Marie Curie Re- 
search Training Network FLAVIAnet under contract No. 
MRTN-CT-2006-035482 



[1] B. Badelek et al. Int. J. Mod. Phys. A 19 (2004) 5097- 
5186 

[2] I. F. Ginzburg, G. L. Kotkin, V. G. Serbo and V. I. Tel- 
nov, Nucl. Instrum. Meth. 205 (1983) 47; 
I. F. Ginzburg, G. L. Kotkin, S. L. Panfil, V. G. Serbo 
and V. I. Telnov, Nucl. Instrum. Meth. A 219 (1984) 5. 

[3] E. Boos et al. Nucl. Instr. Meth. A534 (2004) 250; hep- 
ph/0403113 

[4] A. Pukhov, hep-ph/0412191 

[5] S. Jadach, Z. Was, R. Decker and J. H. Kuhn, Comput. 

Phys. Commun. 76, 361 (1993). 
[6] D. A. Anipko, M. Cannoni, I. F. Ginzburg, A. V. Pak, 



O.Panella, Nucl. Phys. B Proc. Suppl. 126 (2004) 354- 
359; hep-ph/0306138; hep-ph/0410123. 
[7] E. Boos and T. Ohl, Phys. Lett. B 407 (1997) 161, hep- 
ph/9705374. 

[8] I.F. Ginzburg, GL. Kotkin, S.L. Panfil, V.G Serbo, 

Nucl. Phys. B 228 (1983) 285. 
[9] M. Baillargeon, G. Belanger and F. Boudjema, Phys. 

Lett. B 404, 124 (1997)[arXiv:hep-ph/9701368]; Nucl. 

Phys. B 500, 224 (1997) [hep-ph/9701372]; [hep- 

ph/9405359]. 

[10] I.F. Ginzburg, GL. Kotkin. Eur. Phys. J. C 13 (2000) 
295-300 



11 



